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ABSTRACT 

In this paper we investigate the effects that a dynamic dark energy component dominant in the 
universe at late epoc hs has on reionization. We follow the evolution of HII regions with the 
an alytic approach o f Furlanetto & Oh (2005|jn two d ifferent universes for which we assume 
the lPeebles & Ratral (120031) and Brax & Martin (2000) quintessence models and we compare 
our results to the ACDM scenario. We show that, for a fixed ionization efficiency, at the 
same cosmological epoch the morphology of bubbles is dominated by high-mass objects and 
the characteristic size of the ionized regions is slightly smaller than in the ACDM model, 
especially at the latest stages of reionization, due to the higher recombination efficiency. As 
a consequence, the bubbles' 'epoch of overlap' happens earlier than in ACDM. Finally, we 
show how the different evolution of the HII regions affects the transmission of the high-z QSO 
spectra, reducing the Lyman flux absorption at small optical depths. 
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1 INTRODUCTION 

Reionization of the intergalactic medium (IGM) is a crucial epoch 
for the history of the universe, when the neutral gas produced at re- 
combination is ionized by the Ultra- Violet (UV) radiation emitted 
by the first luminous sources. After this stage, the IGM contains 
a small amount of neutral gas, responsible for the absorptions in 
the spectra of background objects. Reionization is still a poorly un- 
derstood process because of the unknown nature of the ionizing 
sources and of the complex physical mechanisms producing the ra- 
diation emission. However, as suggested by recent observations, it 
turns out to be better described as a spatially inhomogen eous and 
not instantaneous phase (see, e.g., ICiardi & Ferran]|200l and ref- 
erences therein). 

While the Lyman-a transmitted flux in the high-z quasar 
(QSO) spectra found by the Sloan Digi tal Sky Survey suggests that 
the end of reionization is at z ~ 6 I Fan et alj|200ll ; iBecker et al.l 
1200 it IWhite et al.l l2003t iFan et alj |2006|) . the last results from 
the Cosmic Microwave Background (CMB) polarization provide a 
Thompson optical depth of r = 0.09 ± 0.03 t hat requires the com - 
pletion of the reionization process at z ~ 10 lSper geletall2007h . 
On the other hand, th e IGM temperature measurements at z < 4 
faui & Haimar]|2003r) show that the reionization epoch occurs at 
6 < z < 10, while the lack of evolution in the Lyman-a galaxy 
luminosity function at 5.7 < z < 6.5 suggests that probably half 
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of the IGM is ionized at z ~ 6.5 dMalhotra & Rhoadsll 20041) . More 
recent estimates, based on high-z HIRES QSOs, seem instead to ar- 
gue against sudd en changes in the IGM properties due to late reion- 
ization at z ~ 6 ( IBecker et al J 20071) : it is worth stressing, however, 
that the use of QSO near zones to pr obe the IGM ionization fr action 
could still be problematic (see, e.g.. lBolton & Haehneltl l2007t) . 

Future observatio ns, such as Lyman-a galaxy surveys 
(Kashikawa et al. 2006), measurements of the CMB polarization 
and of the Sunyaev-Zeldovich effect that will be obtained by the 
Planck satellite, and, more importantly, neutral hydrogen 21 cm 
observations through new generation telescopes (LOFAR, SKA), 
must provide further information to constrain the reionization see- 



cal models (see, e.e.lWyithe & Loebll2003t 1 


Barkana&Loeb 2004; 


Haiman & Holder 2003: Madau et alj|2004 


Choudhury & Ferrara 


2006; iGnedirJ l200d: ICiardi et alj 2003b: 


Wvithe&CerJ 2007; 



Iliev et al.ll2007l) have been developed in order to describe how the 
first sources of UV-radiation impact on the IGM. The basic assump- 
tion of these approaches is to model the relations between the HII 
regions, the ionizing sources that allows to describe the morphol- 
ogy of bubbles, and the galactic physics, such as gas cooling, star 
formation, radiative feedback and radiative transfer. Despite some 
controversial results, the reionization process fits reasonably well 
in a "standard" ACDM cosmology, i.e. in a model where cold dark 
matter (CDM) and baryon density fluctuations grow in a flat uni- 
verse dominated at late epochs by a dark energy component, con- 
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sisting in cosmological constant. The latter is characterized by a 
constant equation of state p = wpc 2 , with w = — 1. 

At present, the best probe of dark energy is provided by SNe 
of type la. Howe ver, even includin g the most recent high— z obser- 
vations (see, e.g. jRiess et alj2007t) . it is still not possible to obtain 
a high-precision determination of the dark energy equation. Tight 
constraints could only be put if strong and unjustified priors on its 
evolution are assumed. While observations suggest that w ~ — 1 
at late epochs, the time evolution of the dark energy component 
is basically not constrained. The so-called quintessence models, 
where w is varying in time, are not excluded. In the CDM cos- 
mologies with dynamic dark energy, the main consequence of hav- 
ing w > —1 at high z is the earlier growth of the matter fluctu- 
ations. This determines a higher number density of haloes in the 
quintessence u niverse than in th e standard cosmology at a fixed 
epoch (see, e.g. lMaio et alj|2006l) . Thereby, the presence of a dark 
energy component might affect the reionization process, requiring 
a different ionization efficiency to fully ionize the IGM at a given 
redshift. 

In this work, we use an analytic approach to investigate how 
reionization proceeds in quintessence cosmologies. In doing this, 
we consider two different cosmological models for which we as- 
sume that the redshift dependence of the equation of state pa- 
rameter w(z) follows the s elf-interacting potentia ls proposed by 
IPeebles & Raid J2003I) and iBrax & Martini feOOCh . The predicted 
scenario obt ained by "painting" th e evolution of the HII regions 
through the jFurlanetto & OrfeOQSt hereafter F05) analytic model, 
is compared to that expected for a ACDM universe. 

The paper is organized as follows. In Section|2]we briefly out- 
line the quintessence cosmologies considered here. In Section|3]we 
review the main features of the F05 model. Section [4] contains the 
results on the evolution of the ionized bubbles and their properties. 
The final Section[5]summarises our main conclusions. 



2 PROBING THE DYNAMIC QUINTESSENCE 
COSMOLOGY 

The main aim of this paper is to investigate the process of cos- 
mic reionization in quintessence cosmologies and compare to the 
predicted scenario for a standard flat ACDM universe. Thus, the 
ACDM cosmology will be our reference case, for which we assume 
that the contributions to the present density parameter from cosmo- 
logical constant, matter, and baryons are Qoa = 0.7, !Jo m = 0.3 
and fiob = 0.046, respectively; the Hubble constant is Ho = 70 
km/s/Mpc (i.e. h = 0.7 in units of 100 km/s/Mpc). We also fix the 
normalization of the power spectrum of the matter fluctuations ac- 
cording to ag — 0.9 and the spectral index is n = 1. These parame- 
ters a re in agreement with the WMAP first-year data jSpergel et al] 
120031) . W e recall that the more recent analysis of the WMAP three- 
year data dSpergel et alj|2007l) suggests slightly different values (in 
particular a lower as and a smaller Oom)- Our parameter choice, 
which is done to allow a direct comparison with the similar analy- 
sis made by F05, can have some small quantitative effect on some 
of the results, but cannot alter the general conclusions of our analy- 
sis, which is aimed at discussing the expected differences between 
models where the dark energy component is provided by a cosmo- 
logical constant or by a dynamic quintessence. 

Thus, the dark energy models we consider are cosmological 
scenarios in which the dynamic dark energy component is charac- 
terized by a self-interacting scalar field $, evolving under the ef- 
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Figure 1. Redshift evolution of the cosmic equation-of-state parameter w 
(top panel) and of the growth factor, given in term of D(z)(l + z) and 
normalized to its value at the present time (bottom panel). Different lines 
refer to ACDM model (solid line), RP model (dashed line) and SUGRA 
model (dotted-dashed line). 



fects of the potential V(&). Here we summ arize the main features 
of these models (see Peebles & Ratra 2003, for more details). 
The potential has to satisfy the Klein-Gordon equation: 



$ + 3H (z)4> + 



. 



(1) 



where H(z) represents the redshift evolution of the Hubble con- 
stant given by the usual Friedmann equation. 

The corresponding ^-evolution of the quintessence parameter 
w is provided by the equation of state 



p _ \ - V(<f) 

P c 2 ii + v($) 



(2) 



such that w — > — 1 if $ 2 <C V(3>). Many analytic expressions have 
been proposed for the potential V($). They describe the present 
dark energy amount simply by setting their amplitude on the initial 
conditions and determine different redshift evolutions for w . In our 
analysis, we use the potentials proposed by Peebl es & Ratral J2003h 
(RP hereafter), 



and bv lBrax & Martini d2000h 



_k_ 

$0! 



y($) = k$- a exp 



<1>- 



(3) 



(4) 



as suggested by supergravity theory (SUGRA hereafter); k has di- 
mension of mass raised to the power (a + 4). The values for k 
and a are fixed by assuming a flat universe with a dark energy 
contribution to the present density parameter flodc = 0.7, and 
Wo = w(z = 0) = -0.85 and -0.83 for SUGRA and RP po- 
tentials, respectively. These choices for wo, even if only marginally 
consistent with the observational constraints (see, e.g. lRiess et al.l 
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Figure 2. The halo distribution at different epochs: redshifts z = 12, z = 9, and z = 6 are shown in the panels, from left to right. For each cosmological 
model, the halo mass function is computed consistently to the PS74 formalism and written in terms of the halo radius. Different curves refer to ACDM model 
(solid line), RP model (dashed line) and SUGRA model (dotted-dashed line). 



l2007t) . have been made with the purpose of emphasizing differ- 
ences with the ACDM model. The remaining parameters have been 
fixed to those of the ACDM model: h = 0.7, erg = 0.9 and n = 1. 
Note that the quintessence models here considered do not violate 
the constraints rec ently obtained from the WMAP three-year data 
dSpergel et all2007t) for the electron scattering optical depth, 

. dL{ z ) 

T e = n e (z)a T — ; — dz , (5) 
Jo dz 

computed at the reionization epoch, assumed to be z r — 6. In the 
previous equation, <tt represents the Thompson scattering cross 
section, and n e (z) and L(z) are the electron density and the co- 
moving distance, respectively. In more detail, we estimate for the 
ACDM, RP and SUGRA models r e = 0.132,0.130 and 0.125, 
which agree (at 1.5cr level) with the WMAP optical depth. Among 
all the possible dark energy models, we concentrate on these two 
since they have been accurately investigated by other authors un- 
der the theoretical point of view and their impact on o bservational 
quant i ties have been extensiv ely addressed (see, e.g.. |Polag et al.l 
120041 : Meneghetti et al. 2005). Furthermore, deviations from the 
ACDM behaviour are larger at high redshifts, in an interesting 
regime for the reionization analysis made here. This is evident 
by looking at the time evolution of the equation-of-state param- 
eter w, shown in the upper panel of Fig. Q] Though the RP and 
SUGRA display similar values for w today, strong differences ap- 
pear at high redshifts. Parametrizing the evolution of w in terms 
of the expansion factor a as w(a) = wo + w a (l — a), we find 
that the RP and SUGRA models can be fitted by w a ~ 0.08 
and 0.55 respectively. These values are still c onsistent with the 
present observational constraints. For example. lLidd!e"e t al. (2006) 
combining data from CMB, SNIa and baryonic acoustic oscilla- 
tions found w a = 0.0 ± 0.8. Of course our results will depend 
on the considered dark energy scenarios: models having parame- 
ters more similar to (different from) the ACDM cosmology (for 
which wo — —l,w a = 0) would show smaller (larger) effects 
on the observables discussed here. These high-redshift differences 
affect the initial phases of the reionisation process. In particular, 
we expect that at the redshifts of interest (z > 6) the RP model 
behaves as an 'intermediate case' between SUGRA and ACDM. 
Note that in quintessence cosmologies, the growth factor D(z) is 
larger than in the standard cosmology, as shown in the lower panel 
of Fig. [T] where the quantity D(z)(l + z) has been normalized 
to its value at the present time. Thus, the main effects of having 
w > — 1 at high redshifts are an earlier formation of the structures 



and a higher abundance of haloes th an in the ACDM cosmology at 
a fixed epoch (see the discussion in lMaio et al . 2006). As an illus- 
trative example, in Fig.[2]we sho w the distribution of the da rk mat- 
ter halo sizes as predicted by the Press & Schechteil d 19741) theory 
(PS74) for the ACDM (solid line), RP (dashed line) and SUGRA 
(dotted-dashed line) cosmologies at three different redshifts. In- 
deed, at the same cosmological epoch, the halo distribution is dom- 
inated by larger objects in the quintessence cases, in particular in 
the SUGRA model. This can strongly affect the ionization process 
of the universe. Indeed nu merical simulations dCiardi et alj|2003l : 
ISokasian et afll2003l. |2004|) showed that reionization is an 'inside- 
out' phenomenon, i.e. it begins in overdense regions and expands 
in underdense regions: this property predominantly originates from 
the fast increase of the abundance of ionising sources. However, 
we remark that a slower (or vanishing) change of the source popu- 
lation should cause a rapid escape of the ionising photons towards 
the external underdense regions: in this case reionization would be 
an outside-in process. 



3 AN ANALYTIC APPROACH TO COSMIC 
REIONIZATION 

To investigate how reionization occurs we use the analytic approach 
proposed by F05. In this section we review the main features of 
the model: for further d etails we refer to the original paper and to 
dFurlanetto et al]2004al hereafter F04). 

3.1 Evolution of bubbles without recombination 

The evolution of the ionized bubbles is determined by the hierarchic 
growth of matter fluctuations. Thi s can be described by making use 
of the extended PS74 formalism dLacev & Colel 1993). A warning 
has to be kept in mind about the use of the Press & Schechter mass 
function, that is proven to not work well for rare hal oes at high 
redsh i fts. In particular, numerical simulati ons (see, e.g.. llliev et al.l 
120061 : iReed et alj|2007l ; iLukic et alj|2007t) show that it underesti- 
mates their abundance by a significant factor, with not negligible 
effects on the bubble distribution. However, we recall that the ex- 
tended formalism has not be developed for mass functions different 
from PS74: for this reason we have to rely on it, even if this can af- 
fect some of the the following results. 

In the PS74 scenario, at a fixed cosmological epoch, an ionized 
bubble grows around a galaxy of mass m ga i Tn m i u (z), where 
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m m i n (,z) represents the virial mass corresponding to T = 10 4 K, 
the temperature at which hydrogen cooling becomes efficient. The 
mass associated to the ionized region is tuhii = C m gai, where Q 
represents the ionization efficiency of the galaxy (here assumed to 
be constant), that depends on the star formation rate, on the escape 
fraction of photons and on the number of HII recombinations. Since 
each region is thought to be isolated, it must contain enough col- 
lapsed mass to fully ionize the inner gas. Thus f co ii(5, m) ^ 
where / co n is the collapsed volume fraction of a region of mass m 
m m i n with an inner overdensity 5. In the F04 formalism, this 
leads to the following condition for the overdensity inside a bubble 
of a given mass m (in Lagrangian space): 



5 m > 5 x {m, z) = 5 c (z) - V2K(C)[a 



2 

min 



,1/2 



(6) 



where K(Q = erf x (l — (" ), a 2 (m) is the variance of density 



fluctuations smoothed on the scale m and a„ 



J (m n 



As 



shown in F04, 5 X represents the ionization threshold for the den- 
sity fluctuations in the Lagrangian space and it is assumed to be 
a linear barrier with respect to a (m): 5 x (m,z) ~ B(m y z) = 
Bo(z) + Bi(z)a 2 (m). Hence it is possible to obtain an analytic 
expression for the distribution of the bubbles with mass in the range 
m ± dm/2: 



n(m, z) = \ - 



2 p I din a I B (z) 



7r m 2 I d In m I a(m) 



exp 



B 2 {m,z) 
2a 2 (m) 



(1) 



where p is the mean com oving matter density of the universe. In a 
similar way, adopting the lLacev & Cold dl993f) formalism, we can 
write the merger rate of the HII regions as: 



>(m 1 ,m T ,t) 
dm^di 



2 lldlnBll dinar 



d In t 1 1 d In mr I 

B(m T , z) 



X 



TV t 1 



exp 



Bl( 



1 



(8) 



where d 2 p(mi , hit , t ) / dm,2dt is the probability per unit time that 
a given halo of mass mi merges with a halo of mass 7712 = mr — 
mi. From equation l[8}, it is possible to define the merger kernel 



Q(mi,m2,t) 



A 2 p{m\,mT,t) 



n(rri2,t) dni2dt 



(9) 



that represents the rate at which each region of mass mi merges 
with a region of mass m,2. Since this quantity suffers from some 
limitations because the asymmetry in its arguments becomes im- 
portant for large masses, then the use of Q S! /m(mi, 7712) = 
l/2[Q(mi, m2) + Q(m2,mi)] is preferred for estimating the 
merger rate of the bubbles. This allows us to define the fractional 
volume accretion for a bubble of mass mi that merges with a mass 
mi: 



V( mi ) — 



v(m 2 ; 

V(mi; 



m,2n(m,2, 



XQsym(mi,ni2, t) 



dtl 
dz\ 



(10) 



Finally, we recall that the global ionized fraction can be calculated 

as Xi = C/coii, g (z), where / co ii,g is the global collapsed volume 
fraction. 



3.2 The recombination-limit effects 

Up to now, the recombination limit has been neglected. As a bubble 
grows, the photons propagate more deeply into the neutral IGM, 
and both the dumpiness and the recombination rate of the ionized 
gas increase. The IGM distribution a nd its ionization state can be 
described using the analytic model of Miralda-Escude et al. 
(MHR00). Analysing numerical simulations at z < 4, they found 
an analytic expression for the volume-weighted distribution of the 
gas density: 



fV(A) = AoA^exp 



(A~ 2/3 - Co) 2 
2(2«5 /3) 2 



(11) 



where A = p/p, 80 is the r.m.s. of density fluctuations smoothed 
on the Jeans mass at fixed z, so So oc (1 + z) , Aq and Co rep- 
resent normalization constants and /3 can be set equal to 2.5 as pre- 
dicted for isothermal spheres. The ionization state of the IGM is de- 
termined by a density threshold A; such that the gas with A < Ai 
is totally ionized and the gas with A > Ai is neutral. Under this 
assumption, the recombination rate can be written as 



A(Ai) = A u / P v (A)A 2 dA = A U C . 



(12) 



where C represents the clumping factor and A u is the recombi- 
nation rate per hydrogen atom in gas at the mean density. In the 
F05 model, A u is assu med consistently with t he A-case of the 
MHR00 model (see also iMiralda-Escudel l2003h : A u oc a A (T), 
where a. 4 (10 4 i<:) = 4x 10" 13 cm 3 s _1 . The MHR00 model also 
provides a relationship between the mean free path Ai of the pho- 
tons and the ionized fraction of gas fV(Ai), namely: 



Ai = A [l--FV(Ai) 



-2/3 



(13) 



where Ao is a normalization constant such that XoH(z) — 60 km 
s" 1 at z < 4. 

In the following, we assume that the mean free path derived in 
the MHR00 model could be used also in the quintessence models, 
since the properties of this function should reflect the gas proper- 
ties at the Jeans scale, which is much smaller than the scales we are 
interested in. However, this approximation deserves further inves- 
tigation with suitable hydrodynamical simulations. 

To consider the recombination process it is necessary to re- 
late the recombination rate to the smoothed matter overdensity. In 
doing this, it must be remarked that the main effect of an inhomo- 
geneous gas distribution is an increasing gas dumpiness and sub- 
sequently an increasing HII recombination rate. As a consequence, 
A u oc (1 + 5). When a bubble grows the ionizing photons are 
able to reach its edge, then the threshold must satisfy the condition 
Ai(Ai) ^ R that sets Ai. However, at the same time, the inner 
high gas dumpiness causes an increase of the recombination rate 
and the photons can be absorbed inside the bubble before reaching 
the edge. Then, for a growing region, the ionization rate has to be 
larger than the recombination rate at every time: 



c df coll (5,R) >AuC{R){1 + s) 



(14) 



where C(R) is computed as in equation l !12t for R = Ai. The 
recombination barrier is obtained by searching for the minimun 5 
in the Lagrangian space that satisfies equation dl4t at each given 
mass. 

The recombination process affects the bubbles geometry. 
When the ionizing photons are totally absorbed by the inner recom- 
bination, the HII region stops growing and reaches a maximum size 
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R max that, in a ACDM universe, slowly increases with decreasing 
redshift, as shown by F05. 

In the excursion-set formalism, the recombination limit has a 
deep impact on the distribution of the bubbles. Assuming that the 
recombination barrier is a vertical line crossing S x at -R ma x, the 
trajectories such that <5(-R ma x) < -B(_R m ax, z) will be incorporated 
into HII regions with m < m max and the mass function reads: 



: (m, z) = 



S(fl max ) 



p{5\R m ^)n(m, z\S, max ■) z)dS . 



(15) 

In the previous equation p(5|_R max ) represents the probability dis- 
tribution at the scale -Rmax for a Gaussian density field 



p(5\R n 



1 



■ exp 



2al, 



(16) 



and n(m, z\5, m. ma x, z) is the conditional mass function for a ran- 



dom walk that begins at (8, <r max ) 



n(m, z\S, m n 



2 p I din a I 
7r m? I d in m I 
q- 2 [-B( m ma x, z) - S] 



exp 



cr 2 ) 3 / 2 

w max y 

[B(m,z) 



2(cr 2 



«) 



(17) 



where <7 max = cr(i? ma x). 
ionization barrier at -R ma x 



Since every trajectory lying above the 
belongs to a bubble with 7? = -Rmax, 
the distribution of such Stromgren regions can be obtained from 
equation l |16t : 



2m n 



-erfc 



B(R n 



(18) 



After the saturation, bubbles can grow only by merging. For a 
single point of the IGM, reionization ends when it is incorporated 
in a recombination-limited region, since the ionizing background 
slightly increases after merging of bubbles. Thus 'overlap' is a lo- 
cal phenomenon. The volume fraction of the IGM in bubbles with 
R > -Rmax is provided by the PS74 formalism using the ionization 
barrier and results to be: 

n(m, z)V (m)dm , (19) 
where V(m) is the volume of the ionized region. 

3.3 The Lyman-Q flux transmission 

The morphology of bubbles affects the absorption of the Lyman- 
q flux emitted from high-redshift sources. Indeed, large ionized 
regions allow the transmission of the Lyman-a forest because 
the extent of neutral regions in the IGM is large enough to re- 
duce the Lyman-Q 'da mping wing' absorption, as pointed out by 
iFurlanetto et al.l j2004ch . As shown by F05, since the recombina- 
tion process affects the late stages of reionization, the way how the 
bubbles saturate can be constrained through the Ly-flux transmis- 
sion. Therefore the observed transmitted flux from high-z galaxies 
can constrain the predicted evolution of HII regions. 

Including the recombination limit in the F05 model allows to 
write the probability of having an optical depth smaller than n for 
the i-th transition, determined by the distribution of the IGM and 



by the morphology of the ionized bubbles, as follows: 



P« n) = 



n(m, z)—dm 

P JO 



Pv(A)dA, (20) 



where mmimin is the minimum mass of the HII regions and A m ax 
is the maximum density for which r < Ti. The analytic expres- 
sion for the inner overdensity of each bubble is obtained by assum- 
ing the equation of state for the polytropic gas T = To A 7 , with 
To = 10 4 K and the ionization equilibrium inside each bubble. Un- 
der these assumptions, the neutral hydrogen fraction can be written 
as a function of the matter overdensity A: 



xhi 



(21) 



where Xe is the correction for the singly-ionized helium and F is 
the ionizing rate per hydrogen atom. It mainly depends on the total 
photons' emissivity er and on the mean free path A, as 



r oc \itT 



V 



3 + V 



(22) 



At the end of reionization er oc £ d/ co n ig /dt, r\ — 3/2 if a starburst 
spectrum is assumed and A; is set to the minimum value between 
the bubble radius and -Rmax - Finally, Pv(A) is thought to be inde- 
pendent of the bubble morphology, that could be a good approxi- 
mation at the end of reionization although we need high-resolution 
simulations to test it. Hence, the relation between the local over- 
density and the IGM optical depth is: 



A(Ti) = 



170 



t] a A (10 4 K) 



Hz) 



C 



d/c, 



3 + r/ a a {To) " v ~ y \^Mpc 

1/(2-0.77) 



A, 



dz 



(23) 



where tgpa is the Gunn & Peterson optical depth for the i-th tran- 
sition. Note that in the equation above the value of the Hubble pa- 
rameter h(z) is taken consistently from the different cosmologi- 
cal models. Finally, the probability for the inhomogeneous IGM to 
have a given optical depth r is obtained by substituting A max in 
equation d20b . 

A warning has to be kept in mind regarding the application of 
our model to derive the Lyman fluxes. Some of the simplifying as- 
sumptions present in the description of the recombination process, 
like the abrupt change of the IG M ionization state as a funct ion of 
its density (see, e.g., discussion in lMiralda-Escude et al.l2000l) . We 
expect, however, that a more realistic treatment would change our 
predictions in the same way for the different cosmological models. 
For this reason, we prefer to present our results in terms of ratios 
between the fluxes derived for the quintessence models and those 
predicted for the ACDM model. 



4 RESULTS AND DISCUSSION 

In this section we present the main results of the application of 
the previous model to cosmological scenarios including a dynamic 
quintessence (see Section|2]l- Notice that F04 and F05 apply their 
model to the 'concordance' ACDM model only. 

First of all, we compute the minimum collapsed mass at each 
cosmological epoch from the mass-temperature relation proposed 
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Figure 3. The ionization barrier is here shown at z = 8, assuming an 
ionization efficiency £ = 6. For each cosmological model, the thick line 
represents S x as defined in equation ^6), while the thin line is its linear ap- 
proximation. The curves refer to the ACDM model (solid line), RP model 
(dashed line) and SUGRA model (dotted-dashed line). 



bv lBarkana & Loebl d200lh : 

where T v i r is the virial temperature of a halo having mass M and /i 
is the mean molecular weight of its inner gas. In order to take into 
account the fact that the IGM inside a HII region is not totally ion- 
ized because of the ionization equilibrium assump tion, here we pre- 
fer to set pi to the mean of the values discussed in Barkana & Loeb 
< l200ll) . Anyway, we checked that fixing the molecular weight to 
the value corresponding to a fully ionized IGM (/1 = 0.6) does 
not significantly affect the model predictions on the observables 
discussed at the end of this section. For example, using fi = 0.6 
changes the Ly-a flux transmission by ~ 1% only, irrespectively 
of the considered cosmological model. 

The values for A c and for the redshift evolution of Q,m(z) are 
computed consistently for the different cosmologies. The presence 
of a dynamic quintessence component affects the hierarchic evolu- 
tion of the ionized regions. Indeed, the earlier growth of matter per- 
turbations causes, at a fixed epoch, the formation of larger ionised 
regions in the quintessence cosmology. Another important conse- 
quence is that, even if the assumption of a linear ionization barrier 
is still correct, this is lower than in the ACDM universe. This is ev- 
ident in Fig. [3] where we show the behaviour of the ionization bar- 
rier in the three different cosmological models, computed at z = 8, 
assuming £ = 6. 

This implies a higher probability that matter fluctuations go 
beyond the ionization threshold. Consequently the density of bub- 
bles increases, resulting in a different morphology of HII regions 
in the three cosmologies considered, as shown in Fig. [4] The earlier 
growth of the matter fluctuations in the quintessence models causes 
a faster evolution of the mass function with respect to the ACDM 
universe and a remarkable increase of the density of the largest re- 
gions. As a consequence, at the same cosmological epoch, we ob- 
tain a higher ionized fraction Xi in the quintessence cases. As an 
illustrative example, we observe that at z = 8.8, reionization is at 
its initial phases in the 'standard' universe (Xi = 0.35) while for 
the SUGRA cosmology this epoch corresponds to its late stages, 
Xi = 0.65. This is almost a factor ~ 2 larger than in the 'standard' 
case. Similarly, at z = 8 reionization is at the final stage in the 
SUGRA universe, Xi = 0.82, while Xi = 0.47 for the 'concor- 
dance' model. We recall that the model we adopt, being based on 



the extended PS74 formalism, can account only in a very approx- 
imate way, most often in the linear limit, for the source cluster- 
ing, which is expected to have important effects on the morphology 
of the HII regions. Since the clustering amplitude depends on the 
source abundances which are different in the models considered 
here, our results could be affected by this bias. Work is in progress 
to properly address this problem by improving our modelling with 
suitable numerical simulations. 

In Fig.|5]we show the different characteristic sizes, i? cha r, of 
the ionized regions, as obtained by fixing ( such that reionization 
ends at 2 = 6 (note that ( is not the same for the three models). 
Neglecting the recombination limit, the effects of the quintessence 
are obvious at the early stages of reionization, since i? c har (repre- 
senting the radius for which the bubbles distribution is maximum) 
is larger in the standard universe than in RP and SUGRA mod- 
els. For example, at Xi ~ 0.2 i? c har = 0.4, 0.6 Mpc for ACDM 
and SUGRA, respectively. But the difference becomes increasingly 
smaller as reionization proceeds. This is caused by the different sat- 
uration regime of the IGM in the quintessence cosmologies. Ionized 
regions are smaller in the RP and SUGRA models at the beginning 
of reionization and their sizes evolve faster than in the ACDM uni- 
verse due to the presence of large neutral voids around them, reach- 
ing the characteristic radius of the standard model in the final stages 
of reionization. 

Since the HII regions grow only by merging after the recom- 
bination limit, it is interesting to investigate how the dynamic dark 
energy component affects the bubbles' merger rates. As an exam- 
ple, Fig. [6] shows the merger probability of a region with mass 
M = 10 14 M o at z = 13 for the ACDM, RP and SUGRA uni- 
verses. Similarly to the ACDM case, also for the quintessence cos- 
mologies the evolution of the bubbles is dominated by merging 
events between large systems, in particular in the late stages of the 
reionization process. The main difference is given by the size of the 
involved regions. For the RP and SUGRA models the merger prob- 
ability is higher with bubbles that are even one order of magnitude 
larger than in the standard universe. 

We can now investigate how the recombination limit due to 
the IGM dumpiness affects the geometry and the evolution of the 
ionized regions in the quintessence universes. As already said, in 
doing this, we assume that the results of the MHR00 simulations 
obtained for a standard cosmology are still valid for a dynamic dark 
energy-dominated universe, since we do not expect large differ- 
ences for the IGM volume-weighted density distribution. As dis- 
cussed above, the dark energy component causes the matter fluc- 
tuations to grow earlier. Hence, since the recombination rate de- 
pends on the inner overdensity of the bubbles, recombination is 
strong already at smaller scales in the quintessence universe, com- 
pared to the ACDM case. Thus, the HII regions reach the equilib- 
rium on scales smaller than in the standard model. This is clear 
in Fig. [7] where we present three different stages of reionization. 
While the ionization threshold does not significantly change be- 
tween ACDM and SUGRA models, the recombination barrier 8 ICC 
extends to smaller scales in the quintessence universes. This effect 
is more prominent at the late stages of reionization process, since 
the bubbles reach the equilibrium on scale of the order of 20 — 30 
Mpc in the ACDM universe instead of the ~ 10 Mpc predicted for 
SUGRA. In dynamic dark energy universes, the 'earlier' (in term 
of comoving scales) recombination barrier involves a smaller value 
for i? m ax, computed as the cross-point of S x and S ICC . The same 
trend was already evident in Fig. [5] where the assumed values for 
( were such that Xi(z = 6) = 1. A peculiarity with respect to the 
standard model is the discontinuous recombination barrier found 
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Figure 4. The morphology of the HII regions. The probability distribution of the bubble sizes is computed here neglecting recombination. In each panel 
(corresponding to different cosmological models), the curves refer to different epochs: z = 8.8 (dotted line), 8.6 (dotted-dashed line), 8.4 (short-dashed line), 
8.2 (long-dashed line) and 8 (solid line). The corresponding values for the ionized volume fraction Xi, computed assuming a ionization efficiency f = 6, are 
also reported. 



100 



10 



0.1 I — ■ — ■ — ■ — ' — ■ — ■ — ■ — ' — ■ — ■ — ■ — ' — ■ — ■ — ■ — ' — ■ — ■ — ■ — I — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — I 
0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 

Xi Xi 
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for the SUGRA cosmology, that crosses the ionization threshold 
more than once. To avoid further complications to the model, we 
choose to set 7? max to the mean value of those achieved at the cross- 
points (which are in any case very close to each other). 

Fig. [8] shows an important effect on the bubble distribution of 
the fact that the maximum radius is smaller, illustrated at z — 8 
for different stages of reionization. As a result, the HII regions with 
M < M max tend to pile up on < 10 Mpc scales in the SUGRA 
universe, in particular at the end of reionization. Furthermore, the 
drop of the ionization threshold causes an increase of the volume 
fraction contained in the recombination-limited regions. 

As mentioned before, for a random point in the IGM, the 
reionization process can be considered complete when the point 
joins a sphere with M — M max . Then, since i? max is smaller in 
the quintessence cosmologies, the density of the points belonging 
to a region with M > M max increases and the volume fraction in- 
side bubbles larger than M max becomes larger. As shown in Fig. [9] 
x ICC increases moving from the ACDM to the SUGRA models as 
reionization proceeds involving different 'epochs of overlapping'. 
In this ~ 0.5 is reached earlier in the quintessence mod- 

els than in the standard universe. This effect is analogous to that 
discussed by F05 for a small mean free path of the ionizing pho- 
tons. 

As illustrated before, the bubble morphology and the IGM ion- 



ization state affect the IGM optical depth r and consequently the 
transmission of the Lyman-a flux. To investigate how the IGM opti- 
cal depth distribution depends on ii max in the quintessence models, 
we compute the probability distribution p(r) of the IGM optical 
depth for the Ly-a, j3, 7 transmission following equation l |20| l. In 
Fig.[l0]we compare p(r) for the RP and SUGRA scenarios to that 
obtained considering the ACDM model, p(t) A cdm, for different 
values of i? max . The transmission is lower for small values of r. 
We find that the trends with i? max for the Lyman-/3 (central panels) 
and Lyman-7 (lower panels) are analogous to that for the Lyman-a 
one (upper panels). 



5 CONCLUSIONS 

The purpose of this work is to give a picture of the reioniza- 
tion epoch in the universes dominated by a dynamic dark energy 
component at late epochs, tracing the HII regions evolution us- 
ing an analytic approach based on the hierarchic growth of mat- 
ter fluctuations. In doing this, we consider two cosmological mod- 
els in which the dark ener gy d ensity varies wi t h tim e driven by 
the lPeebles & Ratral J2003h and iBrax & Martini feOOOh potentials. 
Then, we used the analytic approach proposed by F05 to out- 
line the main differences between the evolution of bubbles in the 
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Figure 6. The hierarchic growth of bubbles. In each panel (corresponding to different cosmological models), we show the fractional volume accretion rate 
of an HII region having mass Mi = 10 14 M Q that merges with a mass corresponding to a given radius R. The results, computed at z = 13, are obtained 
assuming different ionization efficiencies: f = 10 (dotted line), f = 20 (dashed line) and £ = 30 (solid line). 




Figure 7. The recombination limit. In each panel (corresponding to different cosmological models), we show the ionization barrier, computed as in F04 (thin 
curves), and the recombination barrier, computed as in F05 (thick curves). The results are shown at three different stages of reionization: Xi = 0.49 (dotted 
lines), Xi = 0.82 (dashed lines) and Xi = 0.95 (solid lines) at z = 6. 



quintessence models and in the standard ACDM cosmology. Our 
results can be summarized as follows. 

(1) The growth of density fluctuations occurs earlier and the ion- 
ization barrier S x is lower in the RP and SUGRA universes com- 
pared to the ACDM one. This causes a strong increase of the high- 
density regions with respect to the ACDM case at the same epoch. 

(2) Neglecting the recombination limit, the characteristic size of 
the HII regions is smaller in the RP and SUGRA cases at the early 
stage of reionization, but the difference is weakened as reionization 
proceeds. 

(3) In the recombination limit, the early growth of the matter 
fluctuations causes the increase of the IGM dumpiness and the 
inner recombination of the bubbles becomes more efficient. As a 
consequence, the HII regions reach the ionization equilibrium on 
slightly smaller scales and the bubble abundance tends to increase. 



The IGM volume fraction contained in bubbles larger than J? max 
increases requiring an earlier 'epoch of overlap' in the quintessence 
universe compared to ACDM. 

(4) The main effect on the high-z QSO radiation transmission 
due to the different evolution of the HII regions is the lower Lyman 
flux absorption at small optical depths in RP and SUGRA cosmolo- 
gies compared to the ACDM model. 
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Figure 8. The bubble size distribution considering the recombination process. The mass distributions show that in the recombination limit (thick curves) 
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